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By means of first-principles calculations, we studied stable lattice structures and estimated su- 
perconducting transition temperature of CaSi2 at high pressure. Our simulation showed stability of 
the AIB2 structure in a pressure range above 17GPa. In this structure, doubly degenerated optical 
phonon modes, in which the neighboring silicon atoms oscillate alternately in a silicon plane, show 
prominently strong interaction with the conduction electrons. In addition there exists a softened 
optical mode (out-of-plan motion of silicon atoms), whose strength of the electron-phonon interac- 
tion is nearly the same as the above mode. The density of states at the Fermi level in the AIB2 
structure is higher than that in the trigonal structure. These findings and the estimation of the 
transition temperature strongly suggest that higher Tc is expected in the AIB2 structure than the 
trigonal structures which are known so far. 

PACS numbers: 



I. INTRODUCTION 

The calcium di-silicide, CaSi2, has a rhombohedral 
crystal structure at the ambient pressure. In this lay- 
ered structure, superconductivity does not appear down 
to 0.03 K.i In a high pressure phase of the a-ThSi2 type 
structure, which is tetragonal, CaSi2 is known to be- 
come a superconductor with the superconducting tran- 
sition temperature Tc — 1.58 K.^ At the pressure of 
P ~ lOGPa, CaSi2 in the rhombohedral structure un- 
dergoes a structural phase transformation into another 
trigonal structure (Phase III) and superconductivity ap- 
pears with Tc rising up to about 3K4ii In this struc- 
ture, silicon atoms form corrugated honeycomb networks. 
Between two Si honeycomb planes, Ca atoms are inter- 
calated, forming another plane of a triangular lattice. 
Each Ca atom locates just above the center of one cor- 
rugated hexagon formed by six Si atoms in the adjacent 
Si plane. If the pressure increases up to P ~15GPa, an- 
other phase transformation takes place (Phase IV) and 
the corrugated Si planes become nearly flat. This phase 
is called Phase IV. The atomic structure of Phase IV is 
rather close to the AIB2 structure, which has perfectly 
flat silicon planes. Since corrugation remains in Phase 
IV, it is sometimes called the AlB2-like structurei^i^ In 
this structure, Tc further rises up to around 14 K, which 
is the highest record in CaSi2. Here, we note that be- 
fore these experimental findings, the structural transi- 
tion from the trigonal structure to the AIB2 structure 
has been predicted theoretically by one of the authors of 
this study and his coworkers/^ 

The superconductivity in the AIB2 structure attracted 
much attention after finding of a high-temperature super- 
conductor. MgB2r- As for CaSi2, theoretical studies were 
done in low-pressure phases fii^ Satta et al. considered 
possibility of the AIB2 structure in a high pressure con- 
dition, but they could not find this structure with fixed 
cell parameters.^ The electron-phonon interaction and 



the superconducting transition temperature were rarely 
estimated theoretically.^" Thus, CaSi2 has not been stud- 
ied so often in the literature compared with MgB2. This 
another superconductor, however, has its own interest 
and importance. This is because CaSi2 provides us with 
an ideal testing ground on which we can compare several 
polymorphs showing superconductivity. Some of these 
polymorphs resemble each other, but the superconduct- 
ing transition temperature Tc changes its value when the 
structural phase transition takes place. If we could find 
a key factor determining the change of Tc with the struc- 
tural phase transition, it would help us to understand 
this superconducting Zintl-phase compound. 

The purpose of the present study is to clarify the na- 
ture of CaSi2 in high pressure phases above lOGPa. Us- 
ing the first-principles calculations, we did the optimiza- 
tion of the atomic structures in CaSi2 at high pressures. 
We obtained an indication of the pressure-induced phase 
transition from the known trigonal structure to another 
high-pressure phase with the AIB2 structure. Then we 
studied the electronic structures, the phonon dispersion 
relations, and the superconducting transition tempera- 
ture. Those results were compared with those of MgB2. 
We estimated the superconducting transition tempera- 
ture by means of the strong-coupling theory using the 
electron-phonon coupling constants obtained by the first- 
principles calculations. Our results tell us that CaSi2 un- 
dergoes a structural phase transformation to the phase of 
the AIB2 structure and that the new phase is expected to 
be a superconducting phase with much higher transition 
temperature than that of the trigonal structure phase. 



II. CALCULATION METHODS 

In this study, we consider the pressure range of P = 
10 20 GPa. For the structural optimization, we started 
simulations from lattices whose unit cell contain one Ca 



2 



atom and two Si atoms and without any symmetry re- 
quirements. The space-group of the trigonal lattices of 
CaSi2 is PSml and that of AIB2 structure is PG/mmm. 
The WyckofF position of the calcium atom at the la site 
of P3ml is given by (0,0,0), while those of two silicon 
atoms are (1/3,2/3, z) and (2/3,1/3, z) with the inter- 
nal parameter z. When z = 0.5, the structure becomes 
identical to the AIB2 structure. 

For the determination of the electronic structure, we 
utilized the density-functional theorjtiiii^ in the general- 



ized gradient approximation^- with the ultra-soft pseu- 
dopotentials for the atomic potentialJ^ The wave func- 
tions and the electronic charge density are expanded in 
the plane- wave basis. 

For the pseudopotcntials, 3s, 3p and 4s electrons of 
the calcium atom are treated as the valence electrons, 
and for silicon atom 3s and 3p electrons are used. Cal- 
culations are performed by the use of software package, 
the Quantum ESPRESSO 



(a) 



to 6 



optimized — e 
experiment 



(b) 



10 



5 10 15 

Pressure [GPa] 



20 




optimized — e — 
experiment ------- 



(c) 



10 

Pressure [GPa] 



15 



20 



0.5 



0.4 



0.3 



0.2 



0.1 



optimized — e — 
experiment ------- 



5 10 15 

Pressure [GPa] 



20 



FIG. 1: Pressure dependence of lattice parameters: (a) lattice 
constant a, (b) lattice constant c, and (c) internal parameter 
z. The parameters are obtained by the structural optimiza- 
tion at each pressure, and they are indicated by open circles. 
The experimentally observed values^ are represented as closed 
circles. 






FIG. 2: Phonon dispersion of CaSi2. (a) The optimized trig- 
onal structure at 10 GPa. (b) The optimized AIB2 structure 
at 20 GPa. Circles display partial electron-phonon interaction 
A^q, which is explained in the text. 



III. RESULTS OF THE CALCULATIONS 

A. Optimizations of the Structures in High 
Pressure Phases 

The lattice constants and the internal parameter were 
optimized by the constant-pressure variable-cell relax- 



ation using the Parrinello-Rahman method. In this cal- 
culation we used a 12 x 12 x 12 k-point grid in the 
Monkhorst-Pack grid and set the energy cut-off for the 
wave functions at 16 Ry and that for the charge density is 
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FIG. 3: The enthalpy curves of CaSi2 in high pressure. Each 
structure is determined with a fixed value of z. Enthalpy 
values relative to that of the AIB2 structure(z=0.5) is given 
at each pressure. 



at 64 Ry for the trigonal and the AIB2 structure. Though 
these values may be comparatively small, the accuracy is 
enough. This is confirmed by the calculations with larger 
energy cut-offs of 40 Ry and 160 Ry, resulting in almost 
the same optimized structure. 

In Fig. [1] we show the optimized lattice parameters. 
In the pressure range from 10 to 17 GPa, the calculated 
lattice constants are in good agreement with the exper- 
imentally observed values. The relative error of each 
constants, is less than 2%. Above 17 GPa, however, the 
calculated lattice parameters disagree with those of the 
AlB2-like structure. Especially, internal parameter z be- 
comes 0.5 in our calculation while in the experiment^ it 
is less than 0.5 and does not reach that value. Our re- 
sult of optimization shows that CaSi2 becomes the AIB2 
structure above 17 GPa. 

This discrepancy is not due to the pseudopotential 



method adopted in our simulation. We checked the re- 
sults by the use of all electron methods: The full po- 
tential linear muffin tin orbital method which is embod- 
ied in the packaged code developed by S. Y. Savrasov 
and D. Y. Savraso\»ii and the full potential linearized 
augmented-plane wave method which is embodied in the 
WIEN2k code?iS We optimized the structure with con- 
stant cell volume and obtained the same results. In all 
methods, the calculated results indicate that the AIB2 
structure is more stable than trigonal structure at high 
pressure. Here we note that the value of the pressure 
obtained by the first-principles calculation could have an 
error in some cases. In fact, the pressures of pure calcium 
estimated by the generalized gradient approximation cal- 
culation, which is the same method as the present one, 
are much lower than those of experimental values. We 
expect that the AIB2 structure will be observed at higher 
pressures in the experiment. 

To test stability of the AIB2 structure, we calculated 
the phonon frequency in the whole Brillouin zone. The 
density functional perturbation theory was employed^^ 
for the phonon calculation, where 4x4x4 q-point grid 
was used.^^ Phonon dispersion of each phase is shown 
in Fig. [21 We observe that only real frequencies ap- 
pear all over the Brillouin zone, which indicates that the 
AIB2 structure is stable in the pressure range higher than 
17 GPa. 

In Fig. [3l we show pressure dependence of enthalpy 
of some atomic structures. Each structure is given by 
optimizing c/a with fixing the value of z. This graph 
indicates that CaSi2 abruptly transforms from the trigo- 
nal structure (2; = 0.42) to the AIB2 structure(z — 0.5). 
In our simulation, thus no transition from the trigonal 
structure to the AlB2-like structure (z=0.44) occurs. 
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FIG. 4: Electronic band dispersions of the optimized struc- 
tures of CaSi2 in (a) the trigonal structure which is realized 
in the Phase III and (b) AIB2 structure, (c) The electronic 
band structure of MgB2. 



B. The Electronic Band Structure 



tronic band structure of CaSi2 in the AIB2 structures was 



We now analyze the band structures of CaSi2 obtained 
by the Kohn-Sham equations. (Fig. (H^a), (b)) The elec- 
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studied in earlier works. Global features of dispersion 
relations are roughly the same, but we can find some im- 
portant differences between the present results and the 
earlier ones. This is partly due to optimization of the 
crystal structure. In the present work, we have fully op- 
timized structure for the Phase III and the AIB2 phase. 

In Fig. |4] (a), the band structure of the trigonal struc- 
ture is shown. The Si Pz bands, Ca d-band and a p-d hy- 
bridized band touch the Fermi level. Since the Si planes 
are corrugated, Pz bands of Si should rather be called 
"7r(7r*)-like" band. 

In Fig. |4] (b), we show the band structure of CaSi2 in 
the AIB2 structure. First, we can observe d-character in 
some hybridized bands near the Fermi level. One of those 
branches is seen at the F point. We should note that for a 
Ca compound in a high pressure phase Ca d-orbitals often 
appear at the Fermi level. Appearance of the d character 
was pointed out theoretically for CaSi in the CuAu struc- 
ture and CaSis in the CuAua structure.^ Second, the tt* 
band lies also near the Fermi level. At the K point, the 
crossing of the tt and tt* bands is observed, which is a 
characteristic of the AIB2 structure. It looks that the 
electrons are doped in the a* and tt* bands of CaSi2. 
Third, the doubly degenerated p-d hybridized bands at 
the A point are occupied. Along the A-L symmetry line, 
one of those bands becomes almost dispersionless, which 
enhances the density of states around the Fermi level as 
shown below. 

We calculated the pressure dependence of electronic 
density of states at the Fermi level(Fig. [5]). Through 
the structural transition, density of states at the Fermi 
level increases from 0.65 [state/eV] at lOGPa to 1.28 
[state/eV] at 20GPa. As shown in Fig. H in the AIB2 
structure, the s, d and p-d bands go down to Fermi level 
and make electron pockets. This is the reason for the 
enhancement of the density of states. 

Here we compare the electronic band structure of 
CaSi2 with that of MgBz (Fig. g] (c)). Although these 
materials have a close resemblance with each other, the 
band structure of MgB2 has some important differences 
from that of the CaSi2. These AIB2 structures have both 
TT- and cr-bands of sp^ -hybridized orbitals. The a bands 
of MgB2 looks to be partly hole-doped creating small 
two-dimensional hole pockets.^ In CaSi2, on the other 
hand, a bands are fully occupied and a flat a* band lies 
along the A-L line around the Fermi level. The MgB2 
is known to be a two-band system which has tt- and cr- 
bands. CaSi2, however, has additional bands around the 
Fermi level, which are Ca 3d bands. Si 3s band, and p-d 
hybridized bands. 



C. The Electron-Phonon Interaction and the 
Superconductivity 

Assuming the phonon-mediated superconductivity, we 
estimated the superconducting transition temperature 
using the strong coupling theory,^'^'^^ in which the 
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FIG. 5: Pressure dependence of the density of states of CaSi2 
at the Fermi level. Sudden increase happens when the struc- 
ture changes to the AIB2 structure at about 17GPa. 



electron-phonon matrix are calculated by the density 
functional perturbation theory. Our results of the Tc 
shows that when the structural transition occurs, Tc 
rises rapidly and reaches to a value one order of magni- 
tude larger than those in the low pressure phase. The 
estimated superconducting transition temperature are 
shown in Fig. [6l In the trigonal structure, calculated 
results are almost one-tenth of the experimentally ob- 
served values^ which is about 3~4K. This discrepancy 
may be due to utilization of an isotropic approximation 
in the Eliashberg theory. In spite of the low values of the 
estimated T^, we can discuss the pressure dependence of 
Tc. 

Let us examine the origin of the pressure dependence of 
Tc in our theoretical data. According to the McMillan's 
formulated is given by three parameters; the electron- 
phonon coupling constant A, the logarithmic average of 
the phonon frequency ti^iog, and the Coulomb parameter 
fi* , in the following form. 



Tr 



T2 



exp 



1.04(1 + A) 
' X- II* il + 0.62X) 



Here A and Luiog are obtained by the first-principle calcu- 
lations using the density functional perturbation theory. 
As for /i*, we assume the value /i* ~0.1 which holds for 
simple metals. In the present substances, the critical 
temperature is mainly determined by A and wiog. The 
table [J shows A and u!\og for the AIB2 and the trigonal 
structures. Of the two parameters, the more significant 
increase is seen in A. While iviog drops about 10%, the 
increase of the A by a factor of 1.5 in the AIB2 structure 
leads to the enhancement of Tc- 

Here we analyze the electron-phonon interaction. The 
parameter A is given explicitly as follows. 



A = 2 



duj 



a^F{uj) 



using the electron-phonon spectral function, 

^(0) Ek.q \Kl+^?5{0J - ^.q)5(£k)5(ek+q) 



Ekq'^(£k)5(ek+q) 
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Here (0) is the density of electronic states with a single 
spin component at the Fermi level, Wj^q and are phonon 
and electron energies, and M^\^^ is the electron-phonon 
matrix elements. For the mode analysis, we introduce 
partial electron-phonon interaction constant A^q, defined 
by 

^^l^qEkq''5(ek)'5(£k+q') 

from which the mean value is obtained as A = Ej/q ^i^q- 
Using Ai/q, we find the most influential phonon mode 
for the superconductivity and T^, which give large con- 
tribution to the electron-phonon interaction parameter 
A. The contribution is shown in Fig. [2l where Ai,q is 
shown by a circle on each phonon dispersion, and the 
radius is proportional to the contribution to electron- 
phonon interaction parameter A. This figure indicates 
that, in the AIB2 structure, the highest mode at the F 
point is effective. This mode is the mode, in which 
the neighboring silicon atoms oscillate in the anti-phase 
within a Si plane. This feature is the same as that ob- 
served in the MgB2, in which the F,2g mode is the key 
mode of the high-temperature superconductivityj^ We 
can see, however, appearance of the high-frequency peak 
in the phonon density of states and a^F(u;) due to the 
E2g mode of CaSi2 is much similar to the result of AIB2. 



may reduce wiog given by 

wiog = exp J duj — loguj j , 

and does not necessarily work to increase the transi- 
tion temperature as exemplified in iodine. In the case 
of CaSi2 in the AIB2 structure, we observe the high- 
frequency optical branch containing the i?2g mode at the 
F point. In this branch, phonon frequency becomes even 
higher than that of corresponding branch in the trigonal 
structure in the low pressure phase. As a result, wiog is 
kept almost in the same order of magnitude though in the 
low-frequency range the spectral function increases (Fig. 
[2]) . Consequently Tc is not decreased by the phonon soft- 
ening. This means that both Big and phonon modes 
contribute to enhancement of electron-phonon interac- 
tion and Tc- 

TABLE I; Comparison of the electron-phonon interaction pa- 
rameter A and the mean logarithmic frequency ojiog between a 
trigonal structure with the corrugated silicon plane at 10 GPa 
and AIB2 structure with the perfectly flat silicon plane at 
20 GPa. 

structure A a;ic,g [K] 
trigonal 0.27 300 
AIB2 0.41 280 
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FIG. 6: Estimation of pressure dependence of Tc. For our 
simulation CaSi2 takes the structure with the corrugated hon- 
eycomb network of Si atoms in the pressure range from 10 to 
17 GPa, and it transforms to the structure with the perfectly 
flat silicon network above 17 GPa. Sudden enhancement of Tc 
is due to this structural transformation. 



In addition, we see another important mode in the 
AIB2 structure, which is the softened optical mode 
around the F point. This mode is the Big mode corre- 
sponding to out-of-plane motion of silicon atoms, whose 
displacement makes the Si plane corrugated. Due to the 
softening, A,yq has a large value. However, the softening 



IV. SUMMARY 

In this study, we found that the AIB2 structure appears 
as the high-pressure phase of CaSi2, and the enhance- 
ment of the superconducting transition temperature, Tc, 
is expected in the AIB2 phase. The enhancement of Tc is 
due to the enhancement in the electron-phonon interac- 
tion. If we assume that CaSi2 has the phonon-mediated 
superconductivity, the E2g and Big phonon modes play 
an important role in the enhancement of Tc through 
the transformation from the structure with corrugated 
Si plane to the structure with the flat one. 
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